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AN  ANALYTICAL  MODEL  FOR  TURBULENCE-INDUCED 
FLEXURAL  NOISE  IN  LARGE  CONFORMAL  SONAR  ARRAYS 

INTRODUCTION 

Large-area  conformal  arrays  that  can  be  mounted  to  the  hull  of  a  ship 
offer  unique  tactical  advantages  over  towed  arrays.  The  performance  of 
such  arrays  is  usually  limited  by  self -noise  and  platform  noise.  In  the 
latter  category,  the  noise  induced  by  the  turbulent  boundary  layer,  located 
near  the  hull,  is  a  major  concern. 

Boundary  layer  turbulence  produces  a  random  pressure  field  that  is 
detected  by  the  array  as  a  noise  source.  This  is  the  so-called  direct  path 
for  flow  noise.  This  path  for  flow  noise  exists  more  or  less  independently 
of  how  the  sensors  are  supported  and  whether  or  not  they  are  point  sensors 
versus  extended  sensors.  Flow  noise  degrades  the  signal-to-noise  ratio  but 
can  be  reduced  by  using  outer  decoupler  blankets  that  serve  to  attenuate 
the  turbulent  boundary  layer  (TBL)  pressure  field.  The  use  of  extended 
sensors  to  provide  spatial  filtering  of  the  flow  noise  is  also  an 
attractive  way  to  diminish  flow  noise. 

Secondary  sources  for  flow-induced  noise  can  also  be  significant.  If 
the  structural  support  plate  (SSP)  is  relatively  lightweight  and  compliant, 
then  the  TBL  can  induce  flexure  of  the  SSP,  which  then  serves  as  a 
secondary  source  of  noise.  This  noise  can  enter  the  array  via  direct 
flexure  of  the  extended  sensors  or  as  acoustic  noise  radiated  by  the  edges 
of  the  SSP.  The  radiated  component  has  been  addressed  by  other 
investigators  [1-8] .  The  former  source,  flexure  induced  into  the  sensors, 
is  the  focus  of  this  paper. 
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The  problem  will  be  modeled  as  follows:  the  SSP,  sensor  array,  and 
outer  decoupler  (OD)  are  considered  to  constitute  a  curved,  layered 
shell  with  water  on  the  OD  side  and  a  vacuum  on  the  SSP  side.  A  backing  of 
a  vacuum  was  chosen  because  it  is  simple  to  model  and,  in  addition,  it 
represents  a  worst-case  scenario;  i.e.,  the  case  in  which  the  SSP  is  backed 
by  a  pressure  release  baffle.  The  Corcos  [8]  model  will  be  used  for  the 
TBL  pressure  spectrum  although  the  theoretical  development  is  applicable  to 
any  model  of  the  wall  pressure  spectral  density.  This  baseline  model  is 
illustrated  in  Fig.  1.  The  formalism  to  be  developed  makes  no  presumptions 
about  the  boundary  conditions  on  the  plate.  Later  we  shall  assume  that  the 
edges  are  simply  supported  in  order  to  illustrate  a  specific  application. 
Numerous  analytical  studies  of  shell  and  plate  motion  indicate  that,  for  a 
large  shell  or  plate,  simple  supports  usually  give  results  that  accurately 
reflect  the  actual  response  of  a  plate  supported  in  more  complicated  ways. 
The  fluid  loading  on  the  plate  will  be  included  by  using  a  rather  simple 
model  developed  by  Junger  and  Feit  [10] .  The  validity  of  this  model  will 
be  established  by  comparing  the  in-water  displacements  so  derived  with  the 
more  exact  predictions  of  Sandman’s  model  [11]. 
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ARRAY  RESPONSE  TO  PLATB  FLEXURE 

Typically,  the  extended  sensor  array  is  situated  on  or  very  near  the 
SSP  as  shown  in  Fig.  2.  If  the  plate  flexes  in  response  to  an  external 
excitation  the  sensors  will  follow;  therefore,  noise  will  be  generated. 

The  external  excitation  will  also  be  detected  by  the  array  even  if  the 
plate  is  rigid;  this  is  the  so-called  direct  path  for  flow  noise  discussed 
earlier.  In  order  to  analyte  only  the  flexural  contribution,  it  will  be 
assumed  that  the  sensors  do  not  respond  to  the  direct  component. 


Fig.  2  -  Mechanism  by  which  flexural  noise 
is  induced  into  the  array 


The  flexural  response  of  a  piezoelectric  plate  has  been  modeled  by 
Ricketts  [12] .  The  relevant  constitutive  equations  are 


T1  *  °U  S!  *  c?2  S2  -  h31  D3  . 


(la) 
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and 


=  c?«  S,  +  cL  s„  -  h„„  D„  , 

12  1  22  2  32  3  ’ 

(lb) 

t  =  _h31  S1  *  h32  S2  *  ^33  D3  ’ 

(lc) 

T3  =  T4  =  T5  =  0  ’ 

(Id) 

where  T,  S,  E,  and  D  represent  stress,  strain,  electric  field  intensity, 

agd  electric  displacement,  respectively.  The  matrix  components 

C..,  3  ,  and  h. .  are  the  elastic,  dielectric,  and  piezoelectric  material 

lj  33  lj 

constants  [13] .  The  superscripts  indicate  that  the  designated  parameter  is 
held  constant. 

According  to  thin  plate  theory,  the  strains  and  S2  are  related  to 
the  displacement  of  the  plate  as  follows: 


(2) 


where  d  is  the  distance  from  the  sensor  midplane  to  the  neutral  plane. 
and  S2  are  assumed  to  be  constant  through  the  thickness  of  the  sensor 
hydrophone  and  Dg  is  also  constant  through  the  thickness,  as  can  be  shown 
by  applying  Gauss’  law  to  a  dielectric.  Consequently,  Eg  is  constant 
through  the  thickness  and  we  can  write 


(3) 


where  V  is  the  voltage  between  the  electrodes  and  a  is  the  thickness  of  the 
sensor . 


Assume  that  the  hydrophones  are  electrically  connected  in  parallel. 
This  is  equivalent  to  steering  the  array  to  broadside.  Typically, 
piezoelectric  sensors  operate  in  an  open  circuit  mode;  hence,  the  total 
charge  appearing  on  the  electrodes  is  zero.  Therefore,  since  Dg 
corresponds  to  the  charge  density,  we  can  write 
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i  * 

yi 

D3  dx  dy  ,  (4) 

yi 

where  the  integrals  are  performed  over  the  lateral  dimensions  of 
each  sensor.  N  denotes  the  total  number  of  sensors  in  the  array. 

Using  Eq.  (lc) ,  we  obtain 

i  i 

7i 

£+b31  Sl+h32  SJ  **  dy  -  0  •  Cfi) 

yi 


The  first  term  is  independent  of  x  and  y;  therefore, 


Si 


—  dx  dy  = 


Nil 

_ * ..  y . 


(6) 


where  1  and  1  are  the  lateral  dimensions  of  an  individual  sensor, 
x  y 

Combining  Eq.  (5)  with  Eq.  (8)  we  can  write 

V  =  N~T'T  j  I  (h31  Sl+h32  S2]  A  (x’y)  ^  dy  »  <7) 

x  y  -* 

where  A  (x,y)  is  an  array  sensitivity  function.  For  an  array  composed  of 
unshaded  hydrophones  we  can  write 

*  /■_  v\  _  if  (x,y)  is  on  a  sensor 

v  ,y'  lO  if  (x,y)  is  not  on  a  sensor  .  (8) 

Equation  (7)  gives  the  noise  (in  volts)  that  results  from  a  flexural 
response  in  terms  of  the  strain  components  and  Sg.  If  the  external 
excitation  function  is  deterministic,  then  there  will  exist  unique,  well- 
defined  strains  that  can  be  computed  with  thin  plate  or  shell  theories.  On 
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the  other  ha*H,  if  the  excitation  is  a  random  pressure  field  (such  as  TBL) , 
then  the  displacement  and  strains  must  be  thought  of  as  stochastic 
variables,  which  can  be  represented  by  a  probability  distribution  function. 
In  such  cases  the  vu^tage  induced  into  the  array  will  also  be  a  distributed 
variable.  There* jre,  in  order  to  properly  assess  the  noise  due  to  flexure, 
one  must  account  for  the  statistical  nature  of  the  excitation,  which,  in 
our  case,  is  the  turbulent  boundary  layer. 

Before  proceeding  with  the  development  of  a  stochastic  model,  it  is 
convenient  to  express  the  noise  sensed  by  the  array  in  terms  of  am 
equivalent-plane-wave  pressure  field.  This  will  allow  direct  comparison  of 
flexural  noise  with  ambient  sea  noise  and  other  noise  specifications.  In  a 
free  field  environment  the  sensors  operate  in  a  hydrostatic  mode.  An 
incoming  plane  wave  of  amplitude  P  will  produce  an  electric  field,  E^, 
across  the  electrodes  of  a  nonflexing  sensor  of  magnitude 


E3  =  (g33*g31+g32]  P* 


(9) 


where  ggg  is  the  transverse  piezoelectric  constant,  and  g^  and  g^  are  the 
lateral  piezoelectric  constants  for  a  piezoelectric  slab  that  is  poled 
through  its  thickness.  Since  E^  is  constant  through  the  thickness,  the 
voltage  across  the  electrodes  is  V  =  E^a.  Therefore, 


V  = 


*  (^h)1 


(10) 


where  g^  =  g.^  +  g^j  g^2  is  known  as  the  hydrostatic  g  constant. 
Equation  (10)  allows  the  noise,  Eq.  (7),  to  be  expressed  in  terms  of  an 
equivalent-plane-wave  pressure  field  impinging  on  an  array  in  the  free 
field. 
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POWER  SPECTRAL  DENSITY  FOR  FLEIURAL  NOISE 

The  spectral  density  ♦^(w)  for  direct  flow  noise  is  typically  computed 
by  the  following  relationship: 

*pp(")  -  //d2k  A  (£,k*,«)  H  (£,w)  T  (£,*)  P  (k,w)  ,  (11) 

where  A  (i,k  ,u )  is  the  array  function  steered  to  $  ;  H(i,«)  is  the 
hydrophone  function;  T  (k,tf)  is  the  transfer  function;  and  P(k,w)  is 
the  wall  pressure  spectrum.  The  essential  features  of  the  derivation  of 
this  relationship  are  derived  in  Blake  [14]  and  other  references  [15-17] . 
This  equation  gives  the  noise  level  sensed  by  the  thickness  mode  (3-3 
mode);  however,  for  flexure,  the  noise  is  induced  via  an  extensional  mode. 
Therefore,  Eq.  (11)  is  not  appropriate.  Starting  with  the  first  principles 
that  govern  random  vibration  theory,  an  analogous  expression  for  flexural 
noise  can  be  derived.  The  result  will  be  formally  similar  to  Eq.  (11),  but 
the  interpretation  of  the  component  functions  will  be  quite  different. 

The  following  derivation  of  the  power  spectral  density  employs  the 
notation  and  terminology  found  in  Lin  [18] . 

By  combining  Eqs.  (7)  and  (10)  we  can  express  the  equivalent  pressure 
as 

=  Ngh  K  *y  0"  ^31  Sl+h32  S2^  A^X,y^  ***  dy  ' 

As  indicated  previously,  the  excitation  field  is  random;  therefore,  the 
array  output  voltage,  and  hence  the  equivalent-plane-wave  pressure,  must  be 
considered  as  random  variables.  The  transverse  displacement  w  and  the 
lateral  strains  and  S2  are  also  random  variables.  The  power  spectral 
density  for  the  equivalent-plane-wave  pressure  can  be  found  by  taking  the 
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Fourier  transform  of  the  corresponding  correlation  function  Rpp  (t,tf)  that 
is  defined  as 


Rpp(t,t’)  =  E  [P(t)  P(t')]  , 


(13) 


where  E  [  ]  indicates  the  expected  value.  Because  P  depends  linearly 

on  Sj  and  S2,  we  obtain  from  Eq.  (12) 

,2 


Rpp(t»t  )  " 


Vr]  0A£/d2r'[k3iVi(?’t’?',t')  lu) 


+  2  hgi  h32  Rg^s(r,t,r  ,t  )  +  h^  Rg^s(r,t,r  ,t  ) 


A(r)A(r  )  , 


where  &g  3  »  ^g  g  *nd  Eg  g  ,  are  the  cross  correlations  on  strains. 

11  12  2  2 


For  example, 


1  1 


ftg  s  (?,t,?  ,t  )  =  E  [s^ (r , t)  S2 (r*  ,t  )] 
X  2 


(15) 


+  1 

Let  h(r,r  ,t,t  )  denote  the  impulse  response  function  of  the  plate, 
which  is  defined  as  the  displacement  of  the  plate  at  r,t  due  to  an 

impulsive  load  given  by 


Pjir.t)  =  S(t-t  )  <5(t-t  )  ^ 

that  is  applied  at  r  ,t  .  The  principle  of  causality  requires  that 

h(r,  r  ,t,t  )  be  aero  when  t  <  t* .  A  general  excitation  can  be  represented 
by  a  superposition  of  impulses;  likewise,  the  displacement  can  be 
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represented  as  a  superposition  of  impulse  responses  when  the  system  is 
linear.  In  addition,  h  will  depend  only  on  the  difference  t-t' .  Thus,  we 
can  write 

»(r,  t)  =  j  dr  JJ  d2  r  h(r,  r  ,  t-r)  P  (r  ,  r)  ,  (17) 

o  R 

where  P  is  the  applied  pressure  and  R  is  the  region  occupied  by  the  plate. 


The  in-plane  strains  are,  therefore, 


S1(r,  t)  =  -  d  —2=  -  d  J  dr  JJd2  r  h^r,  ?  ,  t-r)P(?  ,  r)  (18) 


o  'R' 


and 


S2  (r,  t)  =  -  d  =  -  d  |  dr  jj  d2  r  h^r,  r  ,  t-r)P(r  ,  r)  , 


(IS) 


where  h^  =  92h/9x2  and  h^  =  32h/8y2. 


Using  Eqs.  (18)  and  (19),  the  correlation  functions  for  strains  [18]  can  be 
written  as 


Vi *• J'’ t<)  ~  d*f0dr  J0dT'  fj i2’  !J i2'  t-r) 


(20) 


*1  +»  •  I  +  + I  I 

hxx  (r  '  8  »  t  ~T  )  Rpp  (*>  T>  8  »  r  )i 


where  R 


PP 


for  Rg 


1S2 


is  the  correlation  on  pressure.  A  similar  expression  is  found 

*ad  Rg  g  .  The  cross  spectral  densities  are  obtained  from  the 
2  2 
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Fourier  transforms  of  the  correlation  functions.  For  example,  the  cross 
spectral  density  for  is  given  by 


q  q  (r» 
S1S1 


■*> 

T 


' }  _  t  T  L  (r’ 


t;  r 


i 

t  ) 


(21) 

„  -i(wt  -  v  t  )  ,  * 

•  exp  v  '  dt  dt 

Assuming  the  excitation  to  be  weakly  stationary,  the  correlation  function 
for  the  wall  pressure  will  depend  only  on  the  temporal  separation 

i 

t-t  .  I**  this  case,  it  can  be  shown  [18]  that  the  correlation 

functions  for  the  responses  also  depend  only  on  the  temporal  separation. 
Subsequently,  the  cross  spectral  densities  depend  only  on  a  single 
frequency  parameter.  Thus,  for  stationary  excitations,  Eqs.  (20)  and 
(21)  may  be  combined  to  yield 


S  (r»  r  »  *0 
S1S1 


d2JJd2  SJR'Jd2  ’  *pp  »>  Hxx<?’  »•  “> 


•  H  ,  ,  (r  ,  s  ,  w)  . 

X  X 


(22) 


Similarly, 


*3^2  -  »>  *  d2/Kf‘)2  *J1'/d2  *  *pp  <*'  »  •  »)  Hxx^’  *• 


(23) 


a  4^ 

•  H  ,  ,  (r  ,  s  ,  *)  , 
7  7 


where  the  frequency  influence  function  H(r,s,v)  is  defined  as 
H  (r,  s,  u)  =  J*  h  (r,  s,  t)  e^^  dt  , 


(24) 
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and  ♦  is  the  Fourier  transform  of  R  with  respect  to  time.  The 
PP  PP  * 

4  4^ 

subscripts  on  H  denote  partial  derivatives  with  respect  to  r  or  r  . 


For  spatially  homogeneous  excitations,  ♦  will  depend  only  on  the 
spatial  separation.  In  this  case,  we  can  replace  the  integrands  of  Eqs. 
(22)  and  (23)  by  their  spatial  Fourier  transforms  to  obtain 

=  d2£  /  V2’  Gxx(?’  u) 


G*.  ,  (r  ,  t,  u)  d2k  , 


(25) 


x  x 


where 


and 


v  &■>  -  (y2  0  *PP  ■“)«i  . 


£.  •)  -  iEj  «) 


i  £«s  ,2  ■* 
e  d  s  . 


(26) 


Similar  expressions  are  obtained  for  Q  and  9C  c  .  The  function 

*2a2  ®la2 

G  (r,(,w)  is  called  the  sensitivity  function.  It  represents 

the  structural  response  at  point  r  when  the  excitation  is  harmonic, 

having  wave  number  £  and  frequency  v. 


For  linear  structures,  Lin  [18]  and  Strawderman  [3]  have  shown  that 
the  sensitivity  function  may  be  written  as  a  superposition  of  normal  modes. 
Using  the  method  of  Lin,  we  obtain 
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exx$’  ?'  »)  “ 


z 

0=1 


2  a 


5x 


0  for  reR 


Sm  (t)  E  (u)  for  reR 


where  fm(r)  =  normal  mode  m  for  the  plate  in  vacuo,  and 


(27) 


(28) 


The  frequency  response  of  mode  m,  HgOO »  is  defined  as  the  modal 
displacement  of  the  panel  when  the  excitation  is  a  unit  harmonic  pressure 
having  wave  numbers  corresponding  to  mode  a.  The  function  S  (£)  is 
commonly  referred  to  as  the  modal  shape  function.  From  Eqs.  (25)  and 
(27),  we  get 

—  A  /1\  »2j*  ^ 


.  .2\r  8  fm(r)  8  fn(r  )  „ 

*S  S  (r,  r  ,  *)  =  d*  y  - a- - H5 H  (w)  H 

&lsl  dx2  '2  “ 

a,n  ox 

•  0  *„  (?.»)  s,A  K  (?)  A 


(26) 


Similarly, 


♦sjSjf1'  *  .  •>  -  d  Z  ~tr  <00  M 

a.n 


By 


0  V  <?.■)  s.  (?)  s;  (?)  A 


(30) 


Eq.  (29)  also  gives  «  by  replacing  0/8x  with  0/8y. 

S2a2 
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Equations  (29)  and  (30)  give  the  cross  spectral  densities  for  strains 

in  terms  of  the  modal  response  of  the  structure  and  the  nave  number 

frequency  spectrum  of  the  excitation.  The  power  spectral  density  for 

flexural  noise  can  now  be  found  by  combining  the  Fourier  transform  of 
Eq.  (14)  with  Eqs.  (29)  and  (30).  Again  invoking  temporal  stationarity, 

the  time  Fourier  transform  of  Eq.  (14)  yields 


V">  ■  [si^y-]2  0  A  1 1 dV  [h3i 


(31) 


+  2  h31h32  *S1S2(r>r  +  h32  *S2S2<r’r  »*0 


A(r)  A(r  ) 


Combining  Eqs.  (29),  (30),  and  (31)  we  obtain 

•\2 


♦  («)  = 

PPV  J 


Ngh  lx  XyJ  m,n 


‘>1  f !  V5'“> S„*  VE>  A 


■n  *2  ?  n  *2 

-mr 


O  4  O  41 
2  0  fm«  9  fn<r  ) 


o 

2  a_'2 


31 


9x  5x 


(32) 


32  „  2 
0y 


dy 


'  2 


A  (r)  A  (r  ) 


If  we  define 


00 


mn 


•n 


*PP  &  *0  sn  tf) 


s;  (i)  d2k 


(33) 
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s2f  a2f. 


a2f  a2f 


I  n  V  J.  W  J.  V  a  v  a 

,  ,2  _ m  n  .  .  n _ m 

and  &  (r  7  r  )  =  o  *  o  *  “31  ^32  To  o 

0,1  32  8x2  Sx  2  31  32  8x  1  8y2 


82f  82f  «  82f  82f 

,  ,  _ m  n  ,2  _ m  n 

*  h31h32  ^2  ^“2  *  h32  6y2  ~^2 


then  Eq.  (32)  can  be  written  as 

,2 


%  00 -  ktd  X  v> a: « v°  D  A 

*  m.n 


“mn^’  r  )  A  (r)  A  (r  )  . 


(34) 


(35) 


Equation  (35)  is  the  central  result  of  this  paper.  It  provides  the 
formal  connection  between  the  wave  number-frequency  spectral  density  of  the 
excitation  field  and  the  power  spectral  density  for  the  equivalent-plane- 
wave  pressure  sensed  by  the  array.  This  relationship  is  the  analog  of  Eq. 
(11),  the  formula  that  is  used  for  direct  path  flow  noise.  The  term  *pin 
Eq.  (35)  is  the  analog  of  the  term  P(£,v)  in  Eq.  (11).  It  accounts  for^he 
spectrum  of  the  excitation  field  through  the  relationship  given  by  Eq.  (33) , 
where  is  identical  to  P(£,v).  The  last  term  in  Eq.  (35)  is  the  analog 
of  the  product  of  the  array  and  hydrophone  functions  in  Eq.  (11) .  This  term 
accounts  for  the  sise  and  location  of  the  sensors  in  the  array.  The  terms 
Hn(t#)  and  H*  (w) ,  called  frequency  response  functions  (FRF),  are  not  found  in 
Eq.  (11) .  These  terms  account  for  the  modal  response  of  the  plate  to  a 
harmonic  excitation;  that  is,  the  response  to  each  modal  component  of 
pressure  in  the  spectrum  of  the  excitation  field.  The  FRF  carries 
information  about  the  material  properties  of  the  plate  as  well  as  the 
effects  of  water  loading  and  intermodal  coupling. 

The  double  sum  over  modes  in  Eq.  (35)  can  be  thought  of  as  representing 
a  type  of  intermodal  coupling.  Thus,  one  could  have  two  types  of  intermodal 
coupling:  (1)  via  the  water  loading  on  the  plate,  and  (2)  via  the  off- 
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diagonal  terms  in  Eq.  (35) .  These  off-diagonal  terms  arise  from  cross¬ 
correlations  of  normal  mode  strains. 


THE  FREQUENCY  RESPONSE  FUNCTION 

The  modal  frequency  response  function  Ha(«)  can  be  found  by  solving  for 
the  steady  state  motion  of  the  system  when  a  unit-amplitude,  harmonic- 
generalised  force  in  mode  m  is  applied.  For  a  rectangular  isotropic  thin 
plate,  simply  supported  and  fluid  loaded  on  one  side,  Lin  [18]  derives  the 
following  frequency  response  function: 


H»w  =  rV  »k  *  kl  1*  -  /»2 

x  y  l  x  yj 


k2  -  k2  -  k2 
a  mx  ay 


,  &  =  (»x,iy  , 


where  L  and  L  are  the  plate  dimensions,  p  is  the  mass  per  unit  area  of 
the  plate,  D  is  the  plate  flexural  rigidity,  k  is  the  acoustic  wave 
number,  and  p ^  is  the  density  of  the  fluid.  The  modal  wave  numbers  are 
given  by 


km  “  L  .  m  -  °»  1»  2,  - 

xx 


•  m  w 

m  =  m  =  0,  1,  2,  - 

y  J 


Equation  (38)  provides  a  good  model  for  the  frequency  response  function 
when  the  mass  and  stiffness  of  the  isotropic  support  plate  is  much  greater 
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than  those  of  the  sensors  and  the  outer  decoupler.  When  the  support  plate 
is  orthotropic,  and/or  comparable  in  mass  and  stiffness  to  the  other 
components,  then  the  entire  structure  must  be  thought  of  as  a  composite 
layered  plate  or  shell.  In  such  cases,  the  frequency  response  function  may 
be  obtained  using  the  Donnel  shell  theory  as  generalized  by  Dong  [10]  to 
include  layered  shells. 

Consider  an  open  shell  of  N  layers  as  shown  in  Fig.  3.  The  boundaries 
of  the  shell  are  defined  as  x=0,  9- 0  and  x=xq,  The  radius  of  the 

shell  is  R.  Each  layer  is  homogenous  and  isotropic.  The  stress-strain 
relationship  for  the  layer  is 


with 


and 


a (k)  1 

X 

Vk) 

= 

.  . 

A(k) 

A11 

A  00 
a12 

0 

A(k) 

a12 

A00 

A11 

0 

0 

0 

A(k) 

a66 

(k)  _ 

lll 

o  1 

1  - 

(k)  _ 

„Ck)  E(k) 

12  ~ 

1  -  ^ 

1 

(k)  _ 

E<k> 

‘66  “ 

1 

6 

X 

e9 

.  7xB  . 

(38) 


(39) 


where  E  and  u  represent  Young’s  modulus  and  Poisson’s  ratio. 
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1st  LAYER 
2nd  LAYER 


REFERENCE 

SURFACE 

k**  LAYER 

n'01  layer 


Fig.  3  -  Composite  shell  model 


Evaluation  of  the  force-resultant  and  moment-resultant  integrals 
requires  piecewise  integration  through  the  thickness.  The  resulting 
equations  of  motion  are 


where  [f]  is  a  differential  operator  matrix  given  by 


F11  *  C11 


C66  a2  a2 

r2  as2  at2  ’ 


(41a) 


r  —  o  —  - 
1 22  “  ©6  0x2 


c22  a2 

r2  a*2 


(41b) 
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r  =  n  si-  +  fn  4.*>n  1  94  °22  a4 

33  11  8x4  e2  l”13  2D“J  ax3  a*2  *  "?  5* 


(41c) 


2  L*  a£_  ?22  a^.  C22  a2 

E  l  12  ax2  E2  as2  * "?  *  '  i?  ’ 


r  =r  tci2*cB6)  a2 

12  J21  B  dx  90 


(41d) 


r  =  r 

1 13  1  31 


Ss  L  n*  a3 

R  9x  U11  ax3 


._66|  a 


11  R2  J  6x5d2  ’ 


(41e) 


r  a  r 

1 23  1 32 


'22  s_  ^22a^_  i  f„.  1  a3 

i2  M  E3  BE3  '  R  lD12*2D6aJ  ^2js  ■ 


where 


(°ij-  Dij’  Dij)  =  X  ‘Ifltwi)'  2  ('t-'k-j . 

k=l 

1  K  -  iJ] . 


(41f) 


n 

p  =  X  4  (wj  • 
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and  p ^  =  Volume  density  of  layer  k,  R  =  Shell  radius,  p  =  Area  density  of 
the  composite  plate  and  are  surface  loads  (pressures) .  All  edges 

are  assumed  to  be  supported  by  shear  diaphragms;  i.e., 

along  x  =  0  and  x  =  xq  ,  y  =  w  =  0  and  N^  =  M^  =  0  ,  (44) 


and 


along  8-0  and  8  =  6q  ,  u  =  w=  0  and  N^  =  M^  =  0  .  (45) 


These  boundary  conditions  are  satisfied  exactly  by  choosing  displacement 
functions  as  follows 


u  =s  u  e  ^*l,t  cos  (As)  sin  (n 8)  , 

▼  =  ▼  q  e  sin  (As)  cos  (n 8)  ,  (46) 

and 

w  =  woe^Wtsin  (As)  sin  (n0)  , 

where  u  is  the  displacement  in  the  x  direction,  y  the  displacement  in  the  8 
direction,  and  w  the  displacement  in  the  s  direction.  Also,  s  is  defined 
as  x/R  and 


A 


m  =  0,  1,  2. . . . 


k  =  0,  1,  2 - 


(47) 
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In  general  n  is  not  an  integer.  The  index  k  is  one  more  than  the 
longitudinal  node  lines  along  the  shell. 

Inserting  Eq.  (46)  into  Eq.  (41) ,  we  obtain 


„  >2  66  2  2 
r\,  =  -  C,„  X  ~  u  *  P*  > 


11 


'11 


R 

C, 


r  —  -  C  -  —  rP  + 

r22  "  C66  X  r2  n  +  p*  • 


r  n  ,4  2  fn  h  1  .2  2  °11  4 

r33  "  D11  A  *  p  1°12'D66J  *  “  *  ~i  ° 


2 

R 


n*  ,2  U11  2 
12  *■?". 


11  2 

1 


-  (C12  *  C66)  . 

1 12  R  An, 


and 


writing 


ri3 

m  £is  x 

R  A 

'  D11 

A  * 

D  * 

l12^. 

r23  = 

e2  •' 

D11 

R3 

3  1  1 
n  +  E 

B11 

=  -ril 

a12  " 

''  ri2 

®22 

=  '  r22 

a13  5 

s  "  ri3 

®33 

=  r33 

a23  s 

=  "  r23  . 

X  n“ 


number  o i 

(48a) 

(48b) 

(48c) 

(48d) 

(48e) 

(48f) 

(49) 


The  resonance  frequencies  of  the  shell  can  be  found  by  solring 
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11 

a12 

a13 

12 

a22 

a23 

=  0  . 

(50) 

13 

a23 

a33 

Model  Verification 

The  natural  frequencies  of  a  closed,  two-layered  cylindrical  shell  were 
computed  from  Eq.  (50).  The  results  are  identical  to  Dong’s  results  [10] 
when  the  reference  surface  is  taken  to  be  the  neutral  fiber. 

The  resonance  frequencies  of  various  open  shells  of  one  layer  of 
isotropic  material  have  been  computed  and  compared  to  the  Donnel-Mushtari 
model  values  reported  by  Leissa  [20].  A  partial  comparison.  Table  1,  shows 
excellent  agreement  with  the  Donnel-Mushtari  model  for  a  single-layered, 
isotropic  open  shell. 
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Table  1  -  Lowest  Resonances  of  a  Single  Layer  Open  Shell  [20] 


xol/mR 

Computed  Resonance  (Hz) 

Leissa  (Hz) 

0.10 

n  =  1/3 

14.27821 

14.2782 

0.25 

2.471323 

2.47132 

1.00 

0.9465441 

0.946544 

4.00 

0.422183 

0.422183 

20.00 

5.1433217B-02 

5 . 14333E-02 

100.00 

2.6873313B-03 

2 . 68829E-03 

0.10 

n  =  1/2 

14.28020 

14.2802 

0.25 

2.472816 

2.47281 

1.00 

0.9308990 

0.930899 

4.00 

0.3814159 

.0381416 

20.00 

3 . 6844123E-02 

3 . 68447E-02 

100.00 

2.3361256E-03 

2 . 32640E-03 

0.10 

n  =  2/3 

14.28297 

14.2830 

0.25 

2.47491 

2.47491 

1.00 

0.9103304 

0.910330 

4.00 

0.3378275 

0.337827 

20.00 

2.7429942E-02 

2.74266E-02 

100.00 

3.7624016E-03 

3.74949E-03 

The  results  given  in  Table  1  show  that  the  model  is  consistent  with 
published  results  that  are  known  to  be  valid  for  thin  plate  or  shell 
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applications.  The  key  question  is  whether  or  not  the  model  provides  a  good 

representation  of  a  typical  support  plate/sensor  array/outer  decoupler  for 

a  practical  sonar  array.  An  illustrative  example,  to  be  discussed  further 

in  a  subsequent  section  of  this  report,  consists  of  a  0.0254m  thick 

cylindrical  steel  shell  overlayed  with  a  0.0508m  thick  elastomer  blanket. 

The  lateral  dimensions  are  L  =  4m  and  L  =  2m.  The  radius  of  curvature  of 

x  y 

the  shell  is  5.2  m.  The  frequency  range  of  interest  is  1  kHz  to  10  kHz. 

The  material  parameters  are 

Elastomer  layer:  E  =  1.0X10**  Pa  v  -  0.49  p  =  1200  kg/m^ 

Steel  support:  E  =  2.1X10^  Pa  v  -  0.3  p  -  7800  kg/m^  . 


Many  resonances  are  found  in  the  frequency  range  of  1  kHz  to  10  kHz  as 

shown  in  Fig.  4.  The  same  computation  has  been  performed  for  a  flat  plate 

having  the  same  size  as  the  structure  of  interest.  Figure  5  shows  a 

comparison  of  the  results  for  the  flat  plate  and  the  shallow  shell.  The 

differences  are  significant  only  for  modes  having  m  <  17  and  m  <5. 

x  y 

These  modes  are  below  1  kHz;  so,  to  simplify  the  model,  the  curvature  of 
the  shell  can  be  neglected.  The  same  conclusions  were  found  to  be  true  for 
a  fiberglass  support  plate. 


Fig.  4  -  Modes  of  a  shallow  shell 
23 
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Circumferential  wavenumber 

Fig  5.  -  The  relative  error,  in  frequency, 
incurred  setting  the  shell  curvature  to  zero 


As  reported  by  Dong  [19] ,  the  shallow  shell  oodel  can  be  converted  to  a 
flat  plate  model  by  making  the  following  transpositions: 

R  -*•  • 

V  0  <51) 

R0q  *  Ly,  the  length  along  curved  edge. 

The  equation  of  motion  Eq.  (40)  now  becomes 

^33  '  w  =  %  ’  (52) 

with 
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and 


tn  x 
x 


a  x 


®x  =  0 ,  1 ,  2 ,  .... 
a  =  0,  1,  2,  . 


(54) 


For  a  single  layer  of  isotropic  material  these  equations  reduce  to  the 
classical  equation  of  motion  for  a  thin  rectangular  plate. 


Validity  of  the  Thin  Plate  Assumption 

A  general  criteria  for  the  range  of  validity  of  a  thin  plate  model  does 
not  exist.  The  critical  frequency  proposed  by  Junger  and  Feit  [10]  for  a 
0.0254-mm  thick  steel  plate  is  f  *  8006  Hz,  and  even  higher  for  a 
fiberglass  plate.  Therefore,  a  thin  plate  approximation  would  seem  to  be 
valid  below  about  8  kHz. 


Water  Loading  and  Driving  Force 

Let  qp  denote  the  forcing  function,  and  let  represent  the  acoustic 
pressure  created  in  the  fluid  by  the  plate  motion.  Equation  (52)  becomes 

r33  w  -  *  Pa-  (55) 

We  assume  that  w  can  be  expressed  as  a  sum  of  in  vacuo  eigenvectors 


=  11 


sin 


a  =1 
x 


=1  xy 


Similarly,  qp  and  Pft  can  be  expressed  as: 


(58) 
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and 


<¥  = 


%m  81D 

*  y 


M si0  (V) 1 


p 


a 


sin 


k  x]  sin  fk  y] 

.  “x  J  l  “y  J 


to  obtain  a  linear  set  of  equations. 


(57) 


(58) 


Modal  Decomposition  of  the  Forcing  Function 

The  forcing  function  used  in  flow  noise  computation  is  usually  a  wall 
pressure  spectrum.  The  amplitude  depends  on  the  particular  model  used 
to  represent  the  spectrum.  A  given  wave  number  component  is  equivalent  to 
a  forcing  function  that  can  be  expressed  as 


qp  =  qQ  sin  (kx)  sin  (kyy) 


(59) 


Following  Timoshenko  and  Woinowsky-Krieger  [21] ,  a  modal  component  of 

q_  ,  Bq.  (57) ,  can  be  obtained  by 
x  y 


Water  Loading  Effect 


(60) 


This  problem  has  been  studied  by  many  authors  [22-26] .  Most  of  the 
available  models  have  been  summarized  by  Leibowitz  [27] .  The  major 
complication  arising  from  water  loading  is  that  the  orthogonal  in  vacuo 
eigenmodes  become  intercoupled  by  the  water  loading.  That  is,  the 
interaction  of  the  fluid  and  plate  is  a  function  of  all  the  modes. 


E.E.  MONTGOMERY 


We  shall  adopt  the  Junger  and  Feit  model  as  described  by  Leibowitz 
[27];  this  model  was  developed  for  symmetric  modes  of  the  plate.  However, 
the  analytical  result  is  valid  for  both  symmetric  and  antisymmetric  modes 
(Leibowits  [27] ,  Table  1)  . 

The  Junger  and  Feit  model  for  the  water  loading  can  be  decomposed  into 
modal  components  in  the  same  way  as  was  done  for  an  externally  applied 
pressure  [Eq.  (60)].  The  result  for  mode  (m,n)  is 


where 


I 

mnpq 


pq 


(61) 


i) 


mnpq 


m+n _ 

-  2l  (-i>p"q  v« 


pq 


(62) 


"f 

2 

Cos  7^Lx  Cos 

2  7,V7xd7r 

\s-iH] 

1/2 

K-^J 

[W 

fr  -  iwm  1  L  L  /4 
I  mnpq  mnpql  x  y 


and 


r 

mnpq 


Re 


pmnpq] 


m 

mnpq 


(63) 


By  examination  of  Eq.  (61)  we  see  that  each  mode  (m,n)  is  coupled  to  all  of 
the  other  modes  via  the  term  I 

mnpq 
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Cross  Coupling 

Leibovitz  [27]  has  shown  that  the  cross  coupling  terms  are  much  smaller 

than  the  self-impedance  components  when  k  L  and  k  L  >3,  which  for  w<<w 

4  x  i  y  c 

is  equivalent  to  k  L  and  k  L  »1  (the  thin  plate  criteria) .  Here  k  is 

®XX  a 

the  acoustic  wave  number  w/c.  This  criteria  is  satisfied  for  the  SSP  over 
the  frequency  band  of  interest.  Moreover,  Sandman  [11]  has  shown  that  when 
moderate  structural  damping  is  included,  the  cross  coupling  terms  are 
negligible.  For  reassurance,  cross  coupling  terms  were  computed 
numerically  and  were  found  to  be  negligible  for  the  problem  of  interest. 


Approximate  Water  Loading 

When  cross  terms  are  negligible,  as  is  the  case  for  the  SSP,  the 
equivalent  modal  pressure  on  the  plate  can  be  expressed  as  [10] 


(64) 


where  p ^  is  the  density  of  water  and  k^  is  the  acoustic  wave  number. 

As  a  final  check  on  this  approximation,  some  typical  modes  were 
computed  and  the  results  compared  to  the  more  exact  results  obtained  by 
numerical  integration  of  Eq.  (62) .  A  representative  comparison  is  shown  in 
Table  2  for  8  kHz,  where  the  agreement  is  expected  to  be  the  worst.  The 
differences  between  the  approximate  values  and  the  exact  values  are  found 
to  lead  to  errors  of  no  more  than  1  dB  in  the  final  result. 
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Table  2  -  Water  Loading:  Comparison  of  Approximate 
and  Exact  Methods 


mngr 
*11  11 

I20  20  20  20 
*21  21  21  21 
I30  30  30  30 


Frequency  =  8000  He 


Exact, 
s/m' 


Approximate 

(Ns/m) 

1.368xl07 
7.887xl07 
j5.222xl07 
j 1 . 253xl07 


(1.357xl07 

(5.208xl07 

(2.160xl07 

(3.281xl07 


-  j  1.666x10®) 

-  j  2. 671x1 07) 
+  j  4.963xl07) 
+  j  1 . 162xl07) 


Flexural  Response  of  the  SSP 

Substituting  Eqs.  (S3),  (56),  (57),  and  (64)  into  Eq.  (55)  leads  to  an 
infinite  set  of  uncoupled  equations 


DU  K*k2  ]  *  2  Pl2*2  °66j 

LX  j) 


k2  k2  -  J 

*X 


*  =  q_  -  P 

ID  O  a  B  BO 

x  y  x  y  x  y 


(65) 


The  total  displacement  can  then  be  found  by  summing  over  modes,  as 


m  m 

*  X  X  \«  sin  (V)  "in  M 

a  =1  m  =1  x  7  1  x  J  K  i 


(66) 
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In  practice,  the  response  needs  to  be  determined  for  only  a  finite  number 
of  modes  as  follows: 


M  N 

max  1  max 


*  X  X  v, sin  (V] sin  M  • 


(67) 


mm  mm 


The  minimum  and  maximum  values  must  be  chosen  to  include  all  the  modes  that 
contribute  to  the  response  of  the  plate  in  the  frequency  range  of  interest. 

In  order  to  assess  reasonable  values  of  M  and  N  ,  a  numerical 

max  max 

analysis  was  performed  on  a  prototypical  SSP  subjected  to  a  harmonic 
excitation.  The  following  material  properties  were  assumed: 


Elastomer  layer: 


E  =  l.OxlO6  Pa  V  =  0.49  p  =  1200  kg/m3 


Steel  support: 


E  =  2-lxlO11  Pa  v  =  0.3  p  =  7900  kg/m3 


The  forcing  function  was  taken  as 


P  =  sin  k  x  eiwt  , 
x  ’ 


(68) 


where  k^  =  25  m~*  and  f  =  =  1500  Hs  . 


The  computed  normal  displacement  is  shown  in  Figs.  6  and  7,  taking  into 
account  the  10,  73,  and  886  modes  that  most  contribute  to  the  radiation 
impedance.  In  this  case,  the  radiation  impedance  has  been  computed 
numerically  using  Sandman  *  a  technique.  It  appears  that  about  73  modes  are 
required  to  model  the  essentials  of  the  flexural  response  of  the  SSP.  The 
modes  that  make  the  largest  contributions  are  those  that  have  both  resonant 
frequencies  and  wave  numbers  that  nearly  match  the  excitation.  It  appears 
possible  to  obtain  an  accurate  description  of  the  flexure  by  keeping  only 
the  modes  having  relative  contributions  higher  than  one  percent. 
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Pig 


Pig.  6.  -  The  normal  displacement  of  the  plate 
'  -  1.065  m  computed  with  varying  numbers  of  modes 


7.  -  The  normal  displacement  of  the  plate  computed 
at  x  =  2.13  m  with  varying  numbers  of  modes 
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Finally,  the  flexural  response  as  computed  with  the  radiation  impedance 
from  Sandman’s  [11]  numerical  technique  was  compared  to  a  similar  result 
using  Junger  and  Feit’s  [10]  approximate  formula.  Figures  8  and  8  compare 
the  results  of  these  two  methods  for  two  representative  sections  through 
the  plate;  one  parallel  to  the  x  axis,  the  other  parallel  to  the  y  axis. 

The  two  methods  give  essentially  the  same  results;  i.e.,  the  differences 
are  less  than  the  resolution  of  the  plots.  These  results  support  the  use 
of  the  simpler  Junger  and  Feit  expression  for  the  purposes  of  modeling  the 
TBL  excitation  of  the  SSP.  The  additional  developments  that  follow  will 
assume  the  Junger  and  Feit  model  is  being  employed;  however,  these  could  be 
generalised  in  a  straightforward  way  to  include  Sandman’s  model  for  the 
water  loading.  Modal  coupling,  should  it  be  important,  could  also  be 
accommodated,  but  the  numerical  computations  would  become  much  more 
difficult. 


Fig.  8  -  The  normal  displacement  at  y  =  1.065  m  computed 
with  the  Junger  and  Feit  radiation  term,  compared 
with  the  more  exact  expression  of  Sandman 
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Fig.  9  -  The  normal  displacement  at  x  =  2.13  m  computed 
with  the  Junger  and  Feit  radiation  terms,  compared 
with  the  more  exact  expression  of  Sandman 


Recall  from  the  structure  of  Eq.  (35)  that  the  double  summation  over 
modes  introduces  an  additional  type  of  modal  coupling  into  the  problem. 
Thus,  one  could  have  two  types  of  modal  coupling:  (1)  ria  the  water 
loading  on  the  plate,  and  (2)  via  the  double  summation  in  Eq.  (35)  which 
arises  from  crosscorrelations  between  modes.  We  have  argued  above  that  the 
first  kind  of  modal  coupling  is  negligible  for  our  problem.  We  shall  see 
that  the  second  kind  is  also  negligible. 


The  Modal  Frequency  Response  Function 

The  frequency  response  function  H  («)  is,  by  definition,  the  complex 

s 

modal  amplitude  w^  when  the  applied  excitation  is  a  unit  modal  force  [18] 
given  by 
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A  »v  W  “v  *7 

Fm  =  FlT  Sin  1 -  sin  ~L - 

x  y  x  y 


The  first  factor  serves  only  as  a  convenient  normalisation  factor.  For 
a  composite  panel,  fluid  loaded  on  one  face,  pressure-released  on  the  other 
face,  me  have  already  shown  the  following  relationship  between  modal  force 
and  modal  displacement  [see  Eq.  (65) ] : 


So,  for  the  excitation  of  Eq.  (89),  a  «  TT~L~  *  Thercforo>  by 

x  y  x  y 


definition  of  FRF,  we  have 


THE  ARRAY  FILTER  FUNCTION 


As  derived  previously,  the  eigenfunctions  for  a  plate  of  dimension  L 


and  Ly  are 
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f 

a 


(r) 


sin 


a  r  x 
x 


sin 


m  r  y 

JL _ 1 

L 


y 


(72) 


where  a  represents  a  pair  of  aodal  indices  (a  ,  a  )  that  take  on  integral 

x  y 

values  from  1  to  ».  Proa  Eq.  (34)  we  have 


a 

an 


(r 


♦ 


+ 


fa 


(73) 


h32  kn 


(r) 


(r  )  • 


We  can  perform  the  integration  over  a  single  sensor  located  in  region  R^  e 
7^,  x'i.  j'i)  obtaining 


■  1 

cos  k  x.  -  cos  k  x. 
n  i  n  i 

X  X 

V 

i 

cos  k  y.  -  cos  k  y. 
n  n  7i 

y  y 

k 

A 

k 

n 

n„ 

X 

L  y 

By  suaaing  over  all  regions  R^  occupied  by  sensors  we  obtain  the  total 
array  filter  function. 
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THE  MODAL  SPECTRAL  DENSITY  OF  THE  EXCITATION 


Before  evaluating  ♦  (tf)  it  will  first  be  necessary  to  evaluate  the 


mn 


form  factors  S  (&),  which  are  the  Fourier  transforms  of  the  mode  shapes; 


1.6, 


Sm&  5  J  J  **  sin 

A 

p 


fk  x]  sin  fk  yl  dx  dy 

l  x  J  l  ay  J 


(75) 


where  A^  denotes  integration  over  the  surface  of  the  composite  SSP.  This 
integration  can  easily  be  done  by  parts,  yielding 


s„(£) 


x 


exp[-i(kxLx-»x*)]}  k^  (l  -  exp[-i[kyLy-»yw]]} 


k2  -k2 
"x  x 


k2  -k2 

*V  7 


(76) 


The  power  spectral  density  of  the  excitation  •  (k,  u)  will 

be  represented  by  the  Corcos  [9]  model 


V5’  *o  - 


P  (0,*) 

a^  c 

'2  kc] 

{'2[ 

(k,-kJ 

1 2*  h  kJ 

I2 1  [4  * 

[a2  kc] 

!2]) 

(77) 


where  k  S  u/ U  and  U  is  the  convection  velocity.  In  addition, 
c  c  c 

P(0,w)  2  a+  (1+7)  p\  Vj  /u  , 

where  p^  is  the  water  density  and  is  the  friction  velocity.  The  other 
terms  a+,  7,  a^,  and  Bg  are  empirically  determined  constants. 


We  shall  first  evaluate  the  diagonal  terms  of  ♦  («) ;  that 


mn 


is,  those  terms  for  which  m  =  n.  In  this  case,  we  can  combine  Eqs.  (33), 
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(76) ,  and  (77)  to  obtain 


♦p  w  -  #p  (-)  #p  («)  ;  og  , 

rmm  rm  rm 

*  r 

where 

n  fl  -  cos  fk  L  -  b  t)!  dk 

[  mx  L  l  x  x  x  j  i  x 

\M  5  2  J-  [■5-‘y*ik->Ja*  R  kJ2]  ’ 


7  k^  -  c°»  [k,  \  -  V)]  dk, 

v°  1 2  J-  T^^rPTMi 

and  9 

*1  °2  kc 

Consider  Eq.  (70) .  This  integral  can  be  written  as 


(78) 


(79) 


(80) 


(81) 


4  (*)  =  Real 

Fa 

x 


[ 


2 

X 

[i  - 

.  .‘F,  -  V) 

dZ 

KJ 

2 

[Z*k-J 

2 

z-kc(l+i«i)] 

[z-kc(l-i«1) 

(82) 


where  the  variable  of  the  integration  has  been  replaced  by  conplex  number 

Z.  This  integral  aajr  be  evaluated  using  the  contour  shown  in  Fig.  10a. 

The  numerator  has  been  replaced  by  the  real  part  of  an  exponential  in  order 

to  avoid  a  pole  where  the  Ca  contour  intersects  the  imaginary  axis.  There 

are  simple  poles  at  ~k  and  ♦  k  ,  and  a  second-order  pole  at 

x  “x 

Z=kc(l  ♦  iSj).  Using  the  theorem  of  residues,  we  find  the  following 
results : 
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The  expression  for  4  (v)  can  be  evaluated  in  a  similar  fashion,  the  only 

Pm 

x 

difference  being  that  the  second  order  pole  will  be  located  on  the 
imaginary  axis.  The  result  is 


(84) 


This  completes  the  evaluation  of  f  (*)  for  the  diagonal  terms  m  -  n. 


The  analogous  expressions  for  the  off-diagonal  terms  (a  t  n)  are 
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=  k  k 
m 

y 


fj 

[.  -i  fk  L 

1  -  e  ly  y 

-v)]| 

ifk  L  -n  t] 
1  -  e  l  y  y  y  J 

lK 

[k2-k2  ] 

l  y  "yJ 

fk2-k2 

ly  n> 

,]lk'  *  l‘2kc)2] 

(87) 


These  expressions  are  similar  to  the  previous  ones  except  that  there 
are  now  four  simple  poles  on  the  real  axis  (see  Fig.  10b),  and  the  off-axis 
poles  are  now  simple  poles.  The  evaluation  is  straightforward  if  one 

writes  the  numerator  as  1  +  eMBx  Dx] T  -eH^x^x  nxr]  -e*  (^x^,x-BxT)  . 

The  first  three  terms  may  be  integrated  using  the  upper  contour  in  Fig. 

10b,  whereas  the  fourth  term  may  be  integrated  using  the  lower  contour. 

The  results  will  not  be  presented  because  they  are  at  least  two  orders  of 
magnitude  smaller  than  the  diagonal  terms  for  the  problem  of  interest. 

This  is  to  be  expected,  for  upon  examination  of  Eq.  (33)  it  is  seen  that 

when  tpp(£,v)  is  a  slowly  varying  function  of  £  (as  compared  to  and  Sq)  , 

then  the  contributions  from  the  orthogonal  form  factors  S  and  S*  tend  to 

®  D 

cancel  upon  integration.  Thus,  the  cross  terms  are  negligible  except  when 

L  and  L  are  so  small  that  S  and  S  vary  as  slowly  as  t  (k,w).  Even  if 
xy  mn  PP 

this  were  not  the  case,  the  off-diagonal  terms  in  the  array  filter  function 

are  small  for  the  problem  of  interest.  A  numerical  study  of  the  off- 

diagonal  terms  indicates  that  they  are  typically  eight  orders  of  magnitude 

smaller  than  the  diagonal  terms.  The  largest  contributions  from  off  the 

diagonal  are  terms  for  which  three  of  the  indices  m  ,  a  ,  n  ,  and  n  are 

x  y  x  y 

equal,  and  the  fourth  index  is  only  slightly  different  from  the  others.  In 
such  cases,  the  off-diagonal  contribution  is  at  least  two  orders  of 
magnitude  smaller  than  the  diagonal  terms. 
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(a) 


Fig.  10  (a)  -  Contour  of  integration  for  Eq.  (82) 

<b> 


Fig.  10  (b)  -  Contour  of  integration  for  Eq.  (86) 

Now  consider  the  product  Ha(w)  H*(«).  The  diagonal  terms  have  peaks  at 
frequencies  where  the  real  part  of  the  FRF  goes  to  zero.  Fig.  11  shows  some 
cases  for  light  damping.  At  the  FRF  resonance  frequencies  the  major 
contribution  to  the  modal  sum  comes  from  only  one  mode.  In  the  case  of 
heavy  damping  the  Q  of  these  resonances  is  greatly  reduced;  thus,  more  than 
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one  mode  may  contribute.  Even  so,  the  number  of  modes  that  need  to  be 
retained  is  not  large  because  only  nearest  neighbors  contribute  at  that 
frequency.  For  example,  as  seen  from  Fig  11,  it  is  evident  that  to  compute 
at  about  1-25  one  neec*s  to  retain  mod-,  that  have  indices 

clustered  about  mode  ®x  =  »y  =  nx  =  ny  =  12.  Exclusion  of  the  other  modes 
greatly  reduces  the  computational  complexity  of  Eq.  (33). 


Fig.  11  -  Frequency  response  function  for 
several  modes  for  which  m  =  n 


A  representative  plot  of  off-diagonal  terms  Ha(#)  H*(w)  is  presented  in 
Fig.  12.  It  is  seen  that  the  off-diagonal  terms  may  have  peaks  at  an  FRF 
resonance  frequency,  but  the  amplitudes  of  the  off-diagonal  terms  are  still 
at  least  one  order  of  magnitude  smaller  than  the  diagonal  terms. 
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Fig.  12  -  Comparison  of  FRF  diagonal  terms 
with  FRF  off-diagonal  terms 


NUMERICAL  PROCEDURES 

Because  the  largest  values  of  the  product  BT(«)  H*(w)  occur  at  the 
resonances  of  the  FRF,  an  upper  bound  on  Q(v)  can  be  found  by  evaluating 
Eq.  (35)  at  those  frequencies  alone.  Only  a  fee  modes  about  each  resonance 
need  to  be  retained  in  the  summation.  Additionally,  the  off-diagonal 
contributions  of  each  term  in  Eq.  (35)  are  at  least  two  orders  of  magnitude 
less  than  the  diagonal  terms.  Therefore,  the  product 

V»)  *p  <“>  |  J  V,  *■*  (88) 

rmn 

is  at  least  six  orders  of  magnitude  less  for  m  t  n,  than  for  m  =  n. 
Therefore,  no  off-diagonal  terms  need  to  be  retained. 

The  following  procedure  was  used  to  evaluate  Eq.  (33) .  First,  the 
resonances  of  the  FRF  in  the  frequency  range  of  interest  were  determined. 
There  are  about  20  such  resonances  for  the  cases  studied.  For  each  of 
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these  resonant  frequencies  ippOO  was  computed  by  summing  over  100  of  the 
nearest  neighbor  diagonal  terms.  Some  results  and  conclusions  follow. 


RESULTS 

As  indicated  previously,  the  spectral  density  ♦  («)  represents  the 

rP 

equivalent-plane-wave  pressure  sensed  by  the  array.  In  this  form  Q(w)  can 
be  directly  compared  to  sea  states  and  also  to  specifications  for  flow 
noise.  Three  specific  cases  have  been  studied  numerically.  All  three 
cases  are  dual  layer  laminations  consisting  of  a  SSP  and  an  elastomer  of 
0.051  m  (the  OD) .  The  stiffness  and  mass  of  the  sensor  are  not  included. 
This  is  done  to  obtain  relative  comparisons  that  do  not  depend  on  the 
specific  properties  of  the  sensor.  Thus,  the  three  cases  differ  only  in 
terms  of  the  properties  of  the  SSP.  Figure  13  shows  the  results  for  a 
0.0254  m  steel  SSP  and  a  0.01S1  a  glass  fiber  plate.  Two  values  of  loss 
tangent  mere  used  for  the  glass  fiber  SSP  to  show  the  influence  of  damping. 
The  smaller  value  is  more  representative  of  practical  materials. 


Fig.  13  -  Spectral  density  for  flexural  noise  induced 
into  support  plates  of  various  compositions 
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CONCLUSIONS 

The  results  indicate  that  a  0.0254  m  steel  SSP  of  light  damping  and  a 
0.191  ■  glass  fiber  SSP  of  heavy  damping  will  both  perform  satisfactorily. 
However,  the  glass  fiber  SSP  with  light  damping  may  exceed  many  flow  noise 
specifications.  It  is  also  seen  that  flexural  noise  rises  as  frequency 
decreases.  Therefore,  flexural  noise  arising  from  TBL  excitations  must  be 
considered  if  the  array  is  to  be  used  at  lower  frequencies. 

The  results  presented  here  do  not  account  for  sensor  material 
properties  or  for  the  effects  of  the  inner  decoupler  (ID) .  This  was  done 
for  two  reasons.  First,  the  properties  of  the  sensor  and  inner  decoupler 
vary  with  array  performance  requirements.  In  addition,  it  was  considered 
prudent  to  first  exercise  the  model  for  a  worst-case  scenario  in  order  to 
determine  the  severity  of  a  potential  noise  source.  The  added  stiffness, 
mass,  and  damping  contributed  by  the  sensor  and  the  ID  will  lower  the 
predicted  noise  levels.  Therefore,  the  results  represent  upper  limits  on 
flexural  noise  for  the  various  types  of  SSP’s.  Inclusion  of  the  sensor  in 
the  analysis  is  a  straightforward  process  if  its  properties  are  known. 
However,  to  include  the  ID,  more  analytical  development  will  be  required. 

One  can  draw  a  number  of  conclusions  about  flexural  noise  from  an 
examination  of  Eqs.  (64),  (67),  and  (73)  without  actually  performing 
numerical  computations.  To  minimise  flexural  noise,  one  should: 

*  Minimise  d,  the  distance  from  the  midplane  of  the  SSP  to  the  mid 
plane  of  the  sensor. 

*  Maximize  gh^l3*  rat*os  °*  ^he  hydrostatic  sensitivity 

to  the  lateral  sensitivities. 

*  Maximize  the  number  of  sensors  in  the  array  and  the  lateral 
dimensions  of  the  SSP. 
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•  Minimise  the  gap  between  sensors.  Ve  -ee  from  Eq-  (67)  that  if  the 
gaps  were  closed  in  at  least  one  dimension,  then  there  would  be 
essentially  no  noise  due  to  flexure  of  the  sensors.  However,  noise 
would  still  be  radiated  into  the  water  by  the  edges  of  the  SSP. 

•  Locate  the  neutral  axis  of  the  flexure  at  the  midplane  of  the  sensor. 
For  example,  “the  piesorubber  or  Polyvinlidene  deflouride  (PVDF)  type 
sensors  could  be  mounted  on  both  sides  of  the  SSP.  Then,  their 
flexural  responses  would  cancel.  However,  other  noise  sources  such 
as  direct  path  flow  noise  and  hull  noise  would  be  enhanced  by  such  a 
design. 

SUMMARY 

Large-area,  hull-mounted  conformal  sonar  arrays  are  evolving  in  the 
direction  of  less  weight  and  reduced  hull  standoff.  These  new  array 
designs  employ  lightweight,  flexible,  planar  sensors  and  lightweight 
support  plates.  Increased  levels  of  flow-induced  flexural  noise  may  be  one 
undesirable  consequence  of  replacing  heavier,  stiffer  components  with 
lightweight,  flexible  components.  Models  exist  for  the  radiated  component 
of  flexural  noise.  In  this  report,  a  model  is  developed  that  accounts  for 
a  heretofore  neglected  flexural  noise  mechanism;  i.e.,  direct  coupling  of 
the  flexure  with  the  lateral  displacement  of  a  planar  sensor.  The 
development  of  the  model  was  presented  for  piezoelectric  sensors  such  as 
rubber-lead  titanate  composites  and  PVDF;  however,  the  development  could  be 
adapted  to  consider  other  sensor  types  such  as  fiber  optic  sensors.  The 
expression  derived  for  the  spectral  density  arising  from  this  model  [Eq. 
(35)]  has  terms  of  which  some  are  analogous  to  the  expression  for  the 
direct  path  spectral  density  [Eq.  (11)] .  On  the  other  hand,  some  terms  in 
Eq.  (35)  are  unique;  these  account  for  the  modal  response  of  the  plate  and 
intermodal  coupling  arising  from  fluid  loading  as  well  as  cross-modal 
correlations . 

Some  perspectives  on  the  model,  and  its  implementation,  were  presented 
for  one  particular  array  configuration  that  approximates  arrays  of 
practical  interest.  The  support  plate  with  sensors  and  inner  decoupler  was 
represented  as  a  thin,  finite,  composite  plate,  water  loaded  on  the  sensor 
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side  and  pressure  relieved  on  the  opposite  face.  The  water  loading  was 
modeled  with  an  approximate  expression  derived  by  Junger  and  Peit  [10] . 

The  validity  of  these  assumptions  and  approximations  was  demonstrated  by 
comparing  the  approximated  resonances  and  displacements  with  values 
obtained  by  more  exact  models. 

An  exact  analytical  expression  for  the  spectral  density  of  flexural 
noise  was  derived  for  the  case  where  the  random  excitation  field  is 
characterized  by  the  Corcos  model.  In  general,  this  expression  involves  a 
quadruple  summation,  each  to  infinity,  over  the  normal  modes  of  the  plate. 
However,  it  was  shown  that  only  a  few  modes  need  to  be  included  for  each 
frequency  of  interest.  Therefore,  the  model  is  easily  tractable  for 
practical  applications . 

Finally,  the  sensitivity  of  the  spectral  density  to  various  array 
parameters  was  presented.  The  parameters  of  influence  include  the 
piezoelectric  constants  of  the  sensors,  the  sensor  standoff  distance,  the 
material  constants  of  the  support  plate,  and  the  size  and  spacing 
characteristics  of  the  array. 

Most  of  the  limiting  assumptions  of  the  model  could  be  addressed 
without  modifying  the  basic  approach.  For  example,  if  the  assumption  of  a 
pressure-release  backing  to  the  support  plate  is  not  valid  for  a  given 
application,  then  one  could  employ  a  boundary  condition  that  specifies  a 
known  impedance  which  might  represent,  for  example,  an  ID  placed  between 
the  support  plate  and  the  hull  of  the  ship.  Another  limiting  assumption  is 
that  the  structure  behaves  as  a  thin  plate.  This  assumption  manifests 
itself  in  the  model  via  the  relationships  between  lateral  strains  and  the 
normal  displacement  [Eq.  (2)].  This  limitation  may  be  removed  by  replacing 
Eq.  (2)  with  analogous  expressions  obtained  from  thick  plate  theory. 

The  results  of  this  report  provide  an  analytical  approach  and  a 
mathematical  expression  for  the  noise  induced  in  an  array  of  extended 
sensors  that  derives  from  the  coupling  of  the  lateral  sensitivity  of  the 
sensor  to  the  flexure  of  the  support  structure. 


R.E.  MONTGOMERY 
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